Simple route to non-Gaussianity in inflation 
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We present a simple way to calculate non-Gaussianity in inflation using fully non-linear equations 
on long wavelengths with stochastic sources to take into account the short-wavelength quantum 
fluctuations. Our formalism includes both scalar metric and matter perturbations, combining them 
into variables which are invariant under changes of time slicing in the long- wavelength limit. We 
illustrate this method with a perturbative calculation in the single-field slow-roll case. We also 
introduce a convenient choice of variables to graphically present the full momentum dependence of 
the three-point correlator. 



Introduction — The standard lore of inflation states 
that the statistical properties of its perturbations corre- 
spond to those of Gaussian random fields. Gaussianity is 
sometimes presented as a robust prediction regardless of 
any realisation in the context of a specific model. Such 
a statement may be adequate for comparing inflation 
with topological defects, since the latter are distinctly 
non-Gaussian, but its validity is related to the use of 
linear theory and is only approximate.^ Since gravity 
is inherently non-linear and the scalar fields responsi- 
ble for inflation may have interacting potentials, some 
non-linearity will always be present and manifest itself 
as non-Gaussianity in the cosmic microwave background 
radiation (CMB). Hence the issue is not whether infla- 
tion is non-Gaussian, but how large the non-Gaussianity 
is. In an era of precision cosmology, it is potentially an 
additional observable to be sought in future CMB mis- 
sions. 

The current observational discriminants between in- 
flationary models are the amplitude of the power spec- 
trum, the spectral index, its running, and the tensor to 
scalar ratio all calculated in linear theory. A priori, 
non-Gaussianity offers another quantity, characteristic of 
a model, which can be used as a discriminant. Differ- 
ent models are expected to exhibit a different amount 
and type of non-linearity and hence carry their own non- 
Gaussian signature. Naively, it is expected that any non- 
linearity will be very small given that linear fluctuations 
are small: an order of magnitude estimate is 

However, given the accuracy of future CMB missions, 
one could try to see whether such small non-linearities 
are observable. For example, the Planck satellite will be 
able to detect them, if /nl ~ C'(l) 

The question of non-Gaussianity requires techniques 



Gaussianity also depends on the definition of the initial state of 
the perturbations [J; we shall assume this to be the standard 
vacuum defined at short wavelengths. 



beyond linear theory. Apart from early attempts to cal- 
culate it 3, the issue attracted some attention recently 
after the promise of high quality CMB data from satel- 
lites like WMAP and Planck. The usual theoretical ap- 
proach has been to expand the Einstein equations to 
second order 0| . An extensive survey of such tech- 
niques (including references) is provided in 0. A tree- 
level action treatment of interacting perturbations was 
given in Q. 

In [lol | we presented a new framework for calcu- 
lating non-Gaussianity. It is valid in a very general 
multiple-field inflation setting, which includes the pos- 
sibility of a non-trivial field metric. In this paper we 
illustrate that general formalism by applying it to the 
simple case of single-field inflation. Explicit calculations 
in the multiple- field case can be fomid in jTll |. 

The formalism consists of a long-wavelength approx- 
imation coupled with a stochastic picture for the gen- 
eration of perturbations on the relevant scales. It de- 
scribes the evolution of variables which are invariant un- 
der changes of time slicing in the long- wavelength limit. 
These variables specify the inhomogeneous system, tak- 
ing into account both matter and scalar metric pertur- 
bations, and are defined without recourse to the concept 
of a linear perturbation. The equations are formally very 
simple, yet they encode the full non-linear evolution on 
large scales. These stochastic equations offer a generali- 
sation of currently existing stochastic approaches to the 
description of inflationary fluctuations |l2l |. They are 
suitable for a perturbative analytic treatment which is 
technically less complicated than perturbing the original 
Einstein equations. They are also well-suited as a ba- 
sis for computer simulations which include the full non- 
linear evolution of the perturbations. 

This paper is organised as follows. We start by sum- 
marising the general formalism of 0,0] , simplified to the 
single-field case. Using this we then analytically calculate 
the three-point correlator of the curvature perturbation 
for single-field slow-roll inflation, which would be zero in 
the Gaussian case. Our result is a completely explicit ex- 
pression, containing the full momentum dependence. We 
introduce a convenient choice of variables to graphically 
present this momentum dependence. We also calculate 
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the effects of the decaying mode. We conclude by dis- 
cussing the results. 

Basic formalism — A defining characteristic of in- 
flation is the behaviour of the comoving Hubble radius 
l/{aH), which shrinks quasi-exponentially. A mode with 
comoving wavenumber k is called super-horizon when 
k < aH , and conversely sub- horizon when k > aH. 
The inflaton is taken to be in a vacuum state, defined 
such that sub-horizon modes approach the Minkowski 
vacuum for k ^ aH. After a mode exits the horizon, 
it is described by a classical probability distribution with 
a variance given by the power spectrum |l3l |. We will 
focus attention on the classical super-horizon regime and 
assume that the dynamics on these scales are adequately 
described by dropping from the equations all terms ex- 
plicitly containing second-order spatial gradients 0| (i.e. 
the long- wavelength approximation).^ Then, focusing on 
scalar modes only, spacetime can be described by the 
metric 



(2) 



with a the local scale factor and N the lapse function. 
The local expansion rate is defined asH = a/{Na), where 
the dot denotes a derivative with respect to t. Here, we 
consider inflation to be driven by a single inflaton <j>. We 
also define two local slow-roll parameters by'^ 



e(t,x) 



77(t,x) 



dV/d(l> 
HU 



(3) 



where V is the potential, H = (j)/N and = SttG — 

The fact that we focus on scalar modes only is an 
approximation, since beyond linear theory, in general, 
scalar, vector, and tensor modes will mix. For complete 
consistency the spatial part of the metric in (0) should be 
multiplied by an /i^ matrix with unit determinant, con- 
taining vector and tensor modes. However, in the long- 
wavelength approximation hp would be non-dynamical, 
as shown for example in jl^. Moreover, at second or- 
der in the perturbations (which gives the leading-order 
non-Gaussianity), the scalar modes are only affected by 
the linear vector and tensor modes, the former of which 
are zero, while the latter are subdominant to the scalar 
modes. Hence it seems a good approximation to neglect 



^ Formally this corresponds with taking only the leading-order 
terms in the gradient expansion. We expect higher-order terms 
to be subdominant on long wavelengths during inflation, but this 
statement has only been rigorously verified at the linear level. A 
calculation to higher order in spatial gradients, or, even better, 
a fu ll proof of convergence of the expansion, would be desirable. 
See mill for more details on the validity of the gradient ex- 
pansion beyond linear theory. 

^ These definitions are equivalent to the more standard ones, 
e = -H/{NH'^) and i) = U/{NHn), if the standard equations 
of motion for H and If are used. However, after introducing 
stochastic source terms (as below) this relation changes subtly. 



vector and tensor modes to obtain a leading-order ex- 
pression for non-Gaussianity. Vector and tensor modes 
will be investigated in future work. 

Since we will be dealing with non-linearities, it is use- 
ful to describe inhomogeneity without resorting to the 
concept of a linear perturbation. In this case, it makes 
sense to consider the spatial gradient of quantities. In 
particular we will work with the following combination 
of spatial gradients: 



Ci(t,x) = 9ilna(t,x) 



a,0(t,x), (4) 



which is invariant under changes of time slicing, up to 
second-order spatial gradients Note that, when 

linearised, is the spatial gradient of the well-known C 
from the literature, the curvature perturbation. 

The full non-linear dynamics of the inhomogeneous 
system are described entirely in terms of C^. As discussed 
in ^3] J one finds a second-order differential equation for 
Q , with a set of constraint equations expressing the gra- 
dients of the coefficients in this equation in terms of Q . 
These equations are exact within the long- wavelength ap- 
proximation. Introducing the velocity 9i, we then rewrite 
the equation of motion as a system of two first-order 
differential equations. Finally we introduce stochastic 
source terms in both to take into account the contin- 
uous ffux of short-wavelength quantum modes crossing 
the horizon into the long-wavelength system. We choose 
the gauge where 



t = lii{aH) ^ NH = (1 - e) 



(5) 



since in this gauge horizon exit of a mode, k ~ aH, occurs 
simultaneously for all spatial points. Then the full non- 
linear system of equations for Q in single-field inflation 
can be written as 



; qua 

m 



(i — di — Si 

+ o^ip 



di — Ji 



(6) 



where the source terms Si and Ji are defined below in 
© and the constraints are 



ft In a = —di InH 



1 -i 



0^ 1 

Oi(P = -; zCi, 

K 1 — e 



2i 



ftH = -—Hi (1-? 



(7) 



1 -e 



Since e, fj, etc. depend on Q via the constraints ((TJ, the 
equation of motion © is clearly very non-linear. Note 
that even though these equations contain slow-roll pa- 
rameters, it is not a slow-roll approximation; the slow-roll 
parameters are just short-hand notation for their respec- 
tive definitions lj3Jl and they have not been assumed to 
be small. On the right-hand side of 10, Si and J are 
stochastic source terms given by 
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r d^k 




(27r)3/2 


-K 1 






(27r)3/2 



>V(fc)ifc,e"'"gii„(/c)a(k)+c.c., 
W(fc)ife,e*" (8) 



Qlin{k) 



1 



Qlm(fc) 



Q;(k) + c.c, 



where c.c. denotes the complex conjugate. W{k) is the 
Fourier transform of an appropriate smoothing window 
function which cuts off modes with wavelengths smaller 
than the Hubble radius; we choose a Gaussian with 
smoothing length R = c/{aH) = ce~*, where c~3~5: 



(9) 



In our gauge R and hence W do not depend on Ci- The 
perturbation quantity Qun is the solution from linear the- 
ory for the Sasaki-Mukhanov variable Q = ~a^/2i(/K. 
It can be computed exactly numerically, or analytically 
within the slow-roll approximation (see e.g. Fi- 
nally, a(k) is a Gaussian complex random number satis- 
fying 



(a(k)a*(k')> =(53(k-k') 



(a(k)a(k')) = 0. (10) 



Since Q 9i are smoothed long-wavelength vari- 
ables, the appropriate initial conditions are that they 
should be zero at early times when all the modes are 
sub-horizon. Hence, 



hm Ci = 0, 

t — > — oo 



Hm 6*, = 0. 

t-^ — QO 



(11) 



The full linear solution Qn^ contains both the growing 
and decaying modes. If we neglect the decaying mode 
(see the end of the next section for remarks about the 
validity of this), it was shown in ^ that in the single- 
field case Qiin = NH{1 + e + Tj)Q\in (actually this result 
is true even beyond linear order). Then the system © 
simplifies considerably, because this means that J^i — 0. 
With the initial conditions specified above, we then must 
have that 0i = at all times. Hence we are left with only 



(12) 



Note that for the fully non-linear case this is still non- 
trivial to solve, since Si has a non-linear dependence on 
Ci, as can be seen from in combination with (TJ. We 
can either deal with it directly numerically, or use an 
approximation method to compute it analytically. The 
numerics are the subject of another paper; here we con- 
tinue with the analytic treatment. 

Analytic approximations — In order to analytically 
solve H12|l we have to apply two approximations: in the 
first place we only consider leading-order terms in a slow- 
roll expansion, and secondly we set up an expansion in 



perturbation orders. Hence only in this section do we 
assume that e and rj are small; up to now the results 
were valid for any value. 

To leading order in slow roll, Qiin is given by (see e.g. 

M) 



\in{k) = ^^_^/ZfcPTH^y^(-_fce-*) 
1 



(13) 



with H^y, the Hankel function of the first kind and of or- 
der 3/2. In the second line we have taken the first terms 
in a series expansion for late times; the first term is the 
growing mode and the second the decaying mode. We 
will at first consider only the growing mode, an assump- 
tion which will be justified at the end of this section. The 
overall unitary factor is irrelevant for the correlators and 
will be omitted. 

Equation (|12|l can now be solved perturbatively. At 
first order all quantities in Si take their homogeneous 
background values. Rewriting a = c/{RH), taking H 
and e to be constant (leading-order slow-roll approxima- 
tion), and switching to i? = ce~* as integration variable, 

g-2, 



the end result is (with C 
d^k 



C«(x)=- 



1 Hn 



2V2J (27r)3/2 /c3/2 ^ 



(14) 



Here the subscript Ti. denotes evaluation at horizon cross- 
ing. The time-dependent part e^*"' ^'^^ actually goes to 
the constant value 1 very quickly (in about 3 c-folds af- 
ter horizon crossing), so that C*'^^ is constant on sufh- 
ciently long super-horizon scales, and independent of the 
smoothing parameter c. Taking R we obtain the 
standard power spectrum: 



(15) 



At second order in the perturbations we expand all 
quantities in Si as follows (C = a, e, etc.): 

C(t,x) =C(°)(t) + C(i)(i,x) =c7(°)+a-2c)Xc>^c)(i), 

(16) 

where we use Q to compute diC. In particular we find 



2i (1 - e)0, + 



1 



4c. 



(17) 



Since perturbing Qua only gives next-to-leading-order 
slow-roll contributions (a consequence of the fact that Q 
is the appropriate quantity to use on short wavelengths, 
not C) and in the gauge Q the window function W is un- 
perturbed, all non-linearity is actually contained in the 
aVe factor in the expression © for Si. The resulting 
equation is 



(2) 



(1) 



(18) 
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with solution 



C(^)(x) 



d^kd^k' 1 Hn H-H' 

8 JJ (2^)3 fc3/2fc/3/2^^ 



fc2 + 



fc2+k-k' , 



Ik + k 



— e'('^+'^')-a(k)a(k') (19) 



-^|l-^eH-''')-a(k)a*(k'))+c.c. 



Again, the time-dependent term e~('^^+'° )-R^/2 y^^-y 
quickly goes to 1, so that C^^-* , like C*-^-* , is constant on suf- 
ficiently long super-horizon scales and independent of c. 

From (|19|l we note that (C^^-*) is indeterminate. To re- 
move this ambiguity and also require that perturbations 
have a zero average, we define C = ^ — (C). Expanding 
^ — Q^'^ and switching over to Fourier space, we find 
the three-point correlator (or rather, the bispectrum) to 
be 



;C(xi)C(x2)C(x3))(ki,k2,k3) (20) 
= i^TrfS' (E.k.) [/(ki, k2) + /(ki, kg) + /(k2, kg)] 

with 




(a) 





J 


sqrt(3) 


ki=0 




'ki\ 


k3=k/2/ 


\k2=k/2 






-1 1 
k2=0 ki=k/2 k3=0 

(b) 



/(k,k') 



2 r 



1 i/^ H^, 
16 fc^/c'^ e-H 



{2in' + VH' 



fc2 + 



x^^it^l+(k^k') 



Ik + k 



/|2 



(21) 



This result is independent of c, so that our choice of 
smoothing scale does not matter. Note that it is valid 
to second order in the perturbations, to leading order in 
slow roll, and on sufficiently long super-horizon scales so 
that C has become constant. 

In the limit <^ ki, k2 (and hence ki — — k2 = k), the 
above expression gives (leaving aside the overall factor of 
(2^)'<53(E,k,)): 



(CCC) 



1 



8 k^kl in 



(2e-H3 +m3) 



2\/\/-n, m2\ 



= -n(fc3)(|C(fc)n(IC(fc3)| 



(22) 



with n ~ n—1 the scalar spectral index. From (|21|l . the 
non-linearity parameter /nl has momentum dependence 
in general. In the limit where one of the momenta is 
much smaller than the other two, /nl ^ n during single- 
field slow-roll inflation. Hence, at least in this limit, 
non-Gaussianity is very small in any such inflation model 
which is compatible with observations. The result 122|) 
agrees exactly with the corresponding limit of Q ■ A sim- 
ilar conclusion was also reached in j5j . 

However, since we have obtained the full momentum 
dependence, it is interesting to go beyond this speciflc 



FIG. 1: (a) The three-point correlator lP?7t . for single- 
field slow-roll inflation, multiplied by kik2k'i/{{ki + k2 + 
fci)/2)^/2[(fi:*/16)(i/''/?^)(2? + ^)]-\ plotted to show its de- 
pendence on the relative size of the three momenta, (b) An 
explanation of the triangular domain used, defined in ll23 



limit. Actually the three-point correlator does not de- 
pend on the three full vectors ki,k2,k3, but only on 
three scalar quantities, which can be taken to be the three 
lengths fci , ^2 , fcg (physically this corresponds to statisti- 
cal isotropy). We can redefine variables to get the overall 
magnitude k = ki + k2 + k^ and two ratios 7 and /3, 



k. 



(3 = ~V3 



ki 



which means that 



(23) 



In addition, because ki + k2 



0, 



tions like Iki -I- k2p = fc? and ki • k; 



one can use rela- 
2r — anu Ki • K.2 = (fc| — fc2 — fef)/-^- 
The domain of 7 and (3 is an equilateral triangle as shown 
in figure 1(b) The vertices of the triangle correspond 



to one of the three momenta being zero (the limit for 
(I22|l ). while the sides correspond to one of the momenta 
being equal to half the total sum (fc/2). From its ver- 
tex to the opposite side kg grows linearly. Plotting the 
three-point correlator in such a way demonstrates its 



5 



symmetry most clearly, as shown in figure |l(a)| The 
three-point correlator in that plot has been multiplied by 
fcffc|fcf/((fcJ-|-A:|-|-fc|)/2)^/^ to remove the overah k mag- 
nitude, and the factor of (K^/16)(i/^/e^)(2e -I- 77) has also 
been omitted (assuming that the momentum dependence 
of H-H and the slow-roll parameters can be neglected). 
There is some dependence on the relative magnitude of 
the momenta, but it is only a variation of order unity for 
the single-field case. 

We finish this section by remarking on the decaying 
mode. Even though it disappears very quickly after hori- 
zon crossing, one might think that including it could im- 
prove the accuracy of the result, since the window func- 
tion does operate near horizon crossing. When taking 
into account the decaying mode, jTi is no longer zero, 
and we have to include the 6i equation. Fortunately we 
can still solve this system analytically. In the limit of 
i? — > we find that C''^'' is unchanged, while C^^-* picks 
up an additional term: {2€t-c' +Vw)k'^ / (k'^ + k''^) ctO^') iii 
(|19|l is replaced by 

(25) 

We see that the additional term is not suppressed by 
slow- roll factors, but by a factor and thus depends 
on the smoothing length. However, once we compute the 
three-point correlator, we find that the additional term 
drops out exactly, because of the relative phase factor i. 
Hence for the final result dropping the decaying modes 
is not an approximation. 

Discussion — We have illustrated a new method for 
calculatiiig non-linearity in inflation, which was intro- 
duced in with an explicit calculation in the single- 
field case. Analytic resuls were derived to second order in 
a perturbative expansion and to leading order in slow roll. 
We derived a general expression for the three-point corre- 
lator of the curvature perturbation in Fourier space (the 



bispectrum). In the limit of one of the momenta being 
much smaller than the other two, the expression reduces 
to a simple result: the scalar spectral index times the 
square of the power spectrum. This agrees with previ- 
ous results in the literature but it was derived in a much 
simpler way. In particular, it agrees exactly with the 
tree-level action calculation of Maldacena Q- We also 
demonstrated that the decaying mode makes no contri- 
bution to the three-point correlator at this order. 

We introduced convenient variables to plot the full mo- 
mentum dependence of the three-point correlator. In ^ 
the three-point correlator is also given in the limit of 
all three momenta being equal, which corresponds with 
the centre of the triangle in figure ^ Although similar in 
magnitude, our expression differs in the exact value there. 
This is probably due to our use of a linear stochastic term 
to emulate the effects of sub- horizon perturbations, which 
does not fully include cubic interactions around horizon 
crossing. As is shown in the terms causing this dis- 
crepancy are subdominant in the multiple-field case. 

For the single-field case non-Gaussianity is suppressed 
by slow-roll factors, and hence unobservable for mod- 
els that satisfy the CMB constraints. However, this is 
not necessarily the case for multiple-field inflation mod- 
els [ll|,il9j, where the non-Gaussianity can be large. The 
main strength of our method, apart from being simple, 
is that it is easily applicable to general multiple-fleld in- 
flation; indeed, the basic equations for the multiple-field 
case were already presented in [2, 10]. The formalism is 
also well-suited for numerical implementation; no slow- 
roll approximation is needed, and the end result is a real- 
space realisation, which contains more information than 
just the A^-point correlator. Results from our numeri- 
cal implementation of the formalism will be presented 
in another paper. We believe that the simplicity of this 
formalism, its general applicability, and its suitability for 
numerical simulations make it very useful for future stud- 
ies of non-Gaussianity from inflation. 
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